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Summary. — We review recent progress in determining the infrared singularity 
structure of on-shell scattering amplitudes in massless gauge theories. We present 
a simple ansatz where soft singularities of any scattering amplitude of massless 
partons, to any loop order, are written as a sum over colour dipoles, governed by 
the cusp anomalous dimension. We explain how this formula was obtained, as the 
simplest solution to a newly-derived set of equations constraining the singularity 
structure to all orders. We emphasize the physical ideas underlying this derivation: 
the factorization of soft and coUinear modes, the special properties of soft gluon 
interactions, and the notion of the cusp anomaly. Finally, we briefly discuss potential 
multi-loop contributions going beyond the sum-over-dipoles formula, which cannot 
be excluded at present. 

PACS 11 . 15 . -q - Gauge field theories. 
PACS 12.38.-t - Quantum chromodynamics. 
PACS 12.38.Cy - Summation of perturbation theory. 



1. The role of infrared singularities 

Understanding soft and coUinear singularities is essential for the application of QCD 
to collider physics. Indeed, cross section calculations beyond tree level involve intricate 
cancellations of such singularities in the sum over final states. A detailed understanding of 
the singularities is therefore a prior condition to making precise predictions. Furthermore, 
knowing the singularities, one can resum the dominant radiative corrections to all orders, 
greatly improving the accuracy of the prediction. 

Beyond their immediate significance to phenomenology, infrared singularities open a 
window into the all-order structure of perturbation theory. They admit a simple, iter- 
ative structure, which is common to all gauge theories. Understanding this structure is 
an important step towards understanding scattering amplitudes in gauge theories in gen- 
eral. As an example, recent progress in studying scattering amplitudes in the maximally 
supersymmetric [M — 4) Yang-Mills theory in the planar limit [1-5], has demonstrated 
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that in that case the iterative structure of the amphtude persists in its finite parts. More- 
over, for the first time a bridge was formed between the weak couphng expansion and 
the strong coupling limit. In these studies the infrared singularity structure had a key 
role. In particular, the cusp anomalous dimension jk{o:s) [6-11], which, as we shall see, 
governs soft singularities in any scattering amplitude, was shown to have an important 
role also in determining the finite parts of the amplitude. Today jk{<^s) is the best 
understood anomalous dimension, at both weak [12] and strong coupling [4,13-15]. As 
shown by this example, there is a strong theoretical incentive to gain full understand- 
ing of the singularity structure of scattering amplitudes. Let us now discuss the more 
pragmatic motivation aiming at precision collider phenomenology. 

Cross-section calculations beyond tree level. - The very fact that gauge-theory am- 
plitudes are plagued by long-distance singularities while the corresponding cross sections 
are finite, makes the determination of these singularities essential. The cancellation of 
infrared singularities in cross sections takes place upon summing over degenerate states, 
as originally shown in QED by Bloch and Nordsieck [16], and later proven in [17,18]. 
Virtual gluons generate singularities in amplitudes owing to the integrations over loop 
momenta, which extend over regions where the gluons are soft or collinear with any of 
the hard partons — this puts some internal propagators on shell, leading to singularities. 
In contrast, real emission diagrams arc finite, but singularities appear upon performing 
phase-space integrations over regions where the emitted partons become soft or collinear 
with other partons. The physical cross section is a sum of these two contributions, 
which can be separately computed in dimensional regularization. Schematically, using 
£) = 4 — 2e dimensions with e < 0, one finds cancellations of the form 



where represents a hard energy scale, e.g. the squared centre-of-mass energy s, while 
m|g^ represents the phase-space limit in the integration of the radiated gluon, which 
depends on the observable considered, e.g. a jet mass. Because of their different origin, 
these singularities render any calculation of scattering cross sections beyond tree level 
highly non-trivial. 

At the one- loop order we have a complete understanding of these singularities. This 
forms the basis for general subtraction algorithms, for example based on a colour dipole 
picture [19], rendering the phase-space integration finite. The possibility to perform such 
local subtraction has been invaluable to practical cross-section calculations. Present day 
collider phenomenology requires computations of multi-leg processes in general kinemat- 
ics, in order to allow for maximal flexibility in the application of kinematic cuts dictated 
by the search strategies and experimental needs. This leads to complicated phase-space 
integrations, which can only be done numerically. Thus, a local subtraction of the sin- 
gularities — which guarantees finite integrals — is an absolute necessity. General sub- 
traction algorithms do not exist yet at the multi-loop level and their development is 
of prime importance to precision computations. The first step in this direction is the 
determination of the singidarity structure of amplitudes, the subject of the present talk. 



(1) 





virtual 



real 



Resummation. ~ Beyond fixed-order cross-section calculations, infrared singularities 
also provide the key to resummation of soft and collinear gluon radiation. Singular contri- 
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buttons cancel between real and virtual corrections, but, as shown schematically in eq. (1), 
a residual logarithm survives. These logarithmically-enhanced corrections (Sudakov log- 
arithms) are the dominant radiative corrections for many cross sections. In particular, 
these corrections are parametrically leading when the relevant scales are far apart — for 
example in eq. (1) when m?^^ <C Q^. Whenever the logarithm becomes as large as the 
inverse power of the coupling, it spoils the converges of the expansion, and corrections 
involving powers of hl{m?^^./Q^) need to be resummed to all orders. The situation 
is complicated by the fact that, due to overlapping soft and collinear divergences, each 
order in perturbation theory can give rise to two logarithms, yielding as ln^(m?,j/Q^). 
In this case, resummation is necessary already for \n.{m?^^/Q^) ~ Because these 

logarithms are all generated by the singularities in the amplitude, which always expo- 
nentiate, higher powers of the logarithms at any order in the coupling can be predicted 
based on the singular terms in the first few orders in the loop expansion. This is a key 
ingredient for resummation. 

The most widely used applications of this picture are parton-shower event generators, 
which implement Sudakov resummation with leading logarithmic accuracy, keeping com- 
plete kinematic information on the generated final state. To achieve better precision one 
typically resorts to a more inclusive approach. Indeed, it has been repeatedly demon- 
strated in a variety of applications, e.g. [20-25], that precise predictions can be obtained 
in kinematic regions that are characterized by a large hierarchy of scales upon performing 
Sudakov resummation, provided one gains suflacient control of subleading logarithms and 
related power corrections. 

The theory of Sudakov resummation is especially well developed in inclusive observ- 
ables, where the hard scattering process involves just two coloured partons [22,23,26-31]. 
Such processes are characterized by a single or a double hierarchy of scales. Examples of 
the first category include deep-inelastic structure functions at large Bjorken x [32], and 
Drell-Yan or Higgs production near partonic threshold, or at small transverse momen- 
tum [22,23]. Examples of the second include event shape distributions [20,33], heavy 
quark production [21,34], and inclusive meson decay spectra [24,25]. The Sudakov factor 
in processes involving two (incoming or outgoing) partons, may be written in the generic 
form 

(2) Sud{m^,N)=expla r —\{l-r)^-^ -l]R{m^,r) 

Jo r Lv ^ . ^ J 

V real virtual 

where Ci — Cp or C'a depending on the colour representation of the hard partons 
(fundamental or adjoint), and the radiator is given by 



(3) C,^^ = -'- 
r r 



rm 7, 2 

^IK (a«(fc')) + 2B (a«(rm2)) - 2V {as{r^m'')) 



These two equations summarize, in a compact way, the form of logarithmically-enhanced 
terms in a typical infrared-safe cross section, to all orders in perturbation theory. This 
simple structure is a consequence of factorization, namely the fact that soft and jet 
subprocesses decouple from the hard interaction and are mutually incoherent. Eq. (2) 
incorporates the cancellation between real and virtual singularities anticipated in (1). 
Note that this equation is written in moment space: only then does the real emission 
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phase-space factorise (see, however, [35]). Eq. (3) describes the structure of the Sudakov 
exponent in terms of a few anomalous dimensions which are functions of the running 
couphng only. This additive structure of the exponent is in one-to-one correspondence 
with the phase-space origin of the various corrections: coUinear logarithms, characterised 
by momenta of order m?r, are controlled by B{as); soft (large-angle) logarithms, charac- 
terised by momenta of order mP'r^, are controlled by 'Diets)] finally, the overlap between 
the soft and the coUinear regions is governed by the cusp anomalous dimension (a^), 
which is a universal quantity, the one and only source of double logarithms. 




Singularities in multi-leg amplitudes. - The application of Sudakov resummation to 
hard processes with several coloured partons is less developed, and it will become more 
important for LHC physics. The starting point to perform such a resummation is the 
analysis of the singularity structure of the corresponding scattering amplitudes (fig. 1) 
at fixed angles [36,37], assuming no strong hierarchy between the various kinematic in- 
variants. A priori, upon considering a multi-leg hard process with general kinematics, 
one may expect a complicated singularity structure, more intricate than the simple ex- 
pressions of eqs. (2) and (3). Yet, the goal remains to understand the singularities to 
any loop order in terms on a small set of anomalous dimensions, which are functions of 
the coupling only. 




Fig. 2. - Gluon 'webs' entering the soft function at 2, 3 and 8 loops respectively 
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A further complication arises in a non-Abelian gauge theory, as soon as the amph- 
tude involves more than three hard partons: soft gluon interactions induce correlations 
between kinematic and colour degree of freedom. Soft gluons still exponentiate, but this 
exponentiation now involves matrices defined in a given colour basis (see below). Re- 
summation formulae such as eqs. (2) and (3) may only hold upon diagonalizing these 
matrices. The size of the matrices depends only of the colour representations of the ex- 
ternal hard partons. A priori, however, at each loop order one would expect new colour 
correlations (as suggested by fig. 2), which would require re-diagonalizing the matrix. 

The theoretical understanding of infrared singularities in multi-leg processes has re- 
cently taken a significant leap forward. The first step was taken in [38], where a two- loop 
calculation of the infrared singularities in multi-leg amplitudes was first performed. The 
conclusions were rather suprising at the time: the colour matrix structure of the soft 
anomalous dimension that controls the singularities at two loops turned out to be iden- 
tical to the one at one loop. The next step was taken very recently in [39] and in parallel 
in [40,41]. These papers explained the findings of [38] and proposed a formula that 
generalizes this result to all loops. According to this formula, no new correlations are 
introduced by multi-loop webs. Instead, the correlations generated by soft gluons are 
always described by a sum over two-body interactions between hard partons, and thus 
the matrix structure at any loop order is the same as at one loop. We shall present 
this formula in the next section. This proposal is based on a set of all-order constraints 
(see sees. 6 and 7 below) that relate the singularity structure in any multi-leg ampli- 
tude to the cusp anomalous dimension ^K{ois)- The derivation of these constraints is 
based on factorization and on the universal properties of soft gluon interactions, which 
are described in sees. 3, 4 and 5. Importantly, in amplitudes with four legs or more the 
sum-over-dipoles formula is still an ansatz: although this formula is consistent with all 
available constraints, there may be additional contributions at three loops or beyond, 
which we can constrain but not exclude at present. This issue is briefly summarized in 
sec. 8. 



2. — The sum-over-dipoles formula 

Consider a scattering amplitude Jvl (j),j//i, Q;g(/i^), e) , involving a fixed number n of 
hard coloured partons carrying momenta pi, i = 1 . . . n, all lightlike, p| = 0, and any 
number of additional non-coloured particles. We assume that ultraviolet renormalization 
has been performed (/i being the renormalization scale) thus all remaining singularities 
are associated with long-distances, and can be regularized working in Z? = 4 — 2e dimen- 
sions, with e < 0. The singularities depend on all the kinematic invariants that can be 
formed out of the hard parton momenta, Pi ■ pj (n(n — l)/2 invariants for an n— parton 
amplitude). We work with general kinematics and assume no special hierarchy between 
these invariants; they must all be large compared to Aq^j-,, and their ratios are regarded 
as numbers of order unity. Momentum conservation is not imposed between the coloured 
partons, allowing for any recoil momentum to be carried by non-coloured particles in 
both the initial and final states. Soft and coUinear factorization properties guarantee 
that all infrared singularities can be absorbed into an overall multiplicative factor Z: one 
writes formally 



(4) M{pi/iJ.,as{n'^),e) = Z {pi/ iiF,ois{ii%),e) n{pi/iJ.,iJ./iJ.F,as{tJ-'^)) , 
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where H is finite and can be taken to be independent of e. Note that in general, the 
factorization scale /zf, at which Z is defined, is distinct from the renormalization scale 
II. For simplicity in the following wc choose hf = Eq. (4) should be understood as 
a matrix multiplication in colour space: H is a vector in some colour basis, accounting 
for the hard scattering process, including any loop corrections involving highly virtual 
gluons. These arc necessarily finite. Z is a matrix in this space, mixing the components 
of the vector H, and accounting for soft and coUinear singularities. According to the 
ansatz of Ref. [39] Z assumes the form(^) 



(5) 

Z {pi/ij,,as{ij'^),e) =exp 



A2 



i^x (a.(A^e)) ^ In ( ^^"y'"'^ ) T, • T, 



where the notation ^^^^^ j-^ indiactcs a sum over all pairs of hard partons forming colour 



(a) „(a) 



dipoles, where each pair is counted twice and and Tj • = • T 

where Tj is a generator (^) in the colour representation of parton i. The overall colour 
charge is conserved, 



(6) 



(a) 



The same ansatz was proposed independently in Ref. [40,41]. 

As originally proposed in [43] (see also [44]), singularities in eq. (5) are generated 
exclusively through integration over the D-dimensional running coupling as(A^, e), which 
obeys the renormalization group equation 



(7) 



djjL 



9 OO 

:/3(e,a«) = -2ea,-^E^»(v) 



n=0 



where 6o = (HC^ — 2n/)/3. It is easy to verify that the solution to this equation for 
small coupling and fixed, negative e, is power suppressed at small scales. 



(8) 



' A^ 



,(M^e)+0(a^) 



which guarantees convergence of integrals ranging from A^ = to some fixed scale fj?, 
as in eq. (5). Non-trivial higher- loop corrections enter in (5) only through higher-order 



(^) Following [42] we keep track of the unitarity phases by writing —pi ■ pj = \pi ■ pj\ e"^"^'^, 
where Xij = 1 if i and j are both initial-state partons or are both final-state partons, and Xij = 

otherwise. 

(^) T'"' should be interpreted as follows: for a final-state quark or an initial-state antiquark: 
t^^; for a final-state antiquark or an initial-state quark: — for a gluon: ifcab- 
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corrections to the anomalous dimensions 7]^ (as) and 7j. (ofg). The former — but not 
the latter — is assumed here to admit Casimir scaling, namely to depend on the colour 
representation of the parton i only through an overall factor given by the corresponding 
quadratic Casimir, 



7if (cks) = 2as/T^ -\-0{a'j.) is known explicitly to three-loop order based on the calculation 
by Moch, Vermaseren and Vogt, [12,45]. Potential contributions of higher-order Casimirs 
at four loops and beyond will be briefly discussed in sec. 8. 

The first term in eq. (5), which is governed by the cusp anomalous dimension 7/^ , rep- 
resents the singularities generated by the interaction of large-angle soft gluons described 
by the S function in fig. 1. This term is written as a sum over colour dipoles formed by 
any pair of hard partons; it correlates the kinematic dependence on the Lorentz invariant 
Pi 'Pj with the corresponding product of colour generators, T^ • Tj . This correlation is pre- 
cisely the one present at the one-loop order. This would imply that no new correlations 
are generated by multi-loop webs such as the ones shown in fig. 2 — a highly non-trivial 
statement, which was not yet tested by direct calculations beyond the two-loop level. 

The second term in eq. (5) represents the interaction of coUinear gluons. It is governed 
by the jet anomalous dimension corresponding to each of the external partons, quarks 
or gluons. These anomalous dimensions are defined in (26) below; they depend not only 
on the colour representation of these partons but also on their spin. Their values are 
known to three-loop order, based on the calculation of the quark and gluon form factors 
in Refs. [46,47]; the coefficients have been conveniently collected in Appendix A of [41]. 

Eq. (5) may well be the exact expression for the singularities of any on-shell scattering 
amplitude in massless gauge theories. As already emphasized, the simplicity of this 
result is striking, especially when compared to the lengthy and complicated expressions 
one typically obtains for multi-leg amplitudes. It is also not what one would naturally 
expect looking at the diagrams of fig. 2. Indeed eq. (5) requires that some remarkable 
cancellations take place in these diagrams. It is therefore very interesting to see how 
eq. (5) emerged out of general considerations. This is the goal of the following sections. 

3. Eikonal approximation and rescaling invariance 

The first key ingredient in deriving the constraints on the singularity structure is 
the universal nature of soft gluon interactions, in particular their independence on the 
absolute momentum scale of the hard parton to which they couple. Let us first ex- 
plain the origin of this property and then analyse its consequences in the context of the 
factorization of an on-shcU amplitude. 



Fig. 3. - Gluon bremsstrahlung off an outgoing quark. The final-state quark is on shell, = m^. 

Consider, for example, soft gluon radiation off a hard quark, as shown in fig. 3. We 
assume that after this emission the quark is on shell, so that p^ = (where can 



(9) 



7k ("s) = C^jxias) ; 



Ci = Tj • Tj . 
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vanish, but this is not necessary for the argument that follows). Applying the ordinary 
Feynman rules one obviously obtains a result that depends on the radiating particle spin 
and momentum — the first expression in cq. (10). Considering instead the limit where 
the gluon is soft {k 0) one obtains a much simpler result, 

(10) .(p) (-i,.T( V) , y+^+ri- ^ ^iP)9sT^^^^—, 

\ / (p + fc)^ — + 1£ — » p - k + ie 

which depends only on the colour charge and direction /? of the quark momentum 

(11) g,rw f , . =g.TW- ^ 



p - k + ie f3 ■ k + ie 

Equivalently, we observe that the resulting "eikonal" Feynman rules are invariant with 
respect to rescaling of the quark velocity, (3 — > /t/3, a symmetry property that will be 
central to our discussion in what follows. Finally, we note that the eikonal approximation 
is conveniently formulated by replacing the dynamical hard partons, which provide the 
source for the radiation, with Wilson lines along their classical trajectories. 



(12) $^(0,-oo) = Pexp 



J — cx 



dXp-A{Xp) 



Here rescaling invariance is inherent: it is realised through reparametrization of the 
integral along the path. 

4. — Factorization 

The second key element is the factorization of soft and coUinear singularities in the 
amplitude, illustrated in fig. 1. Following [38,48-51], we write 



/(2ft_^i)2 



(13) X n 



1 — 272— '"«(/^ 



Here the hard function Hk and the amplitude M. l are vectors in the space of available 
color configurations; the soft function Slk is a matrix in this space, while the jet functions 
Ji and Ji do not carry any colour index. The soft matrix S and the jet functions ,/ and 
J contain all soft and coUinear singularities of the amplitude, while the hard functions 
Hk can be taken to be independent of e. Each of the functions appearing in cq. (13) is 
separately gauge invariant and admits operator definitions that are given below. 

The soft and jet functions involve semi-infinite Wilson lines, defined in (12). The 
'partonic jet' function (for, say, an outgoing quark with momentum p) is defined by 

(14) u{p) J (^^^P^,as{ix% = (p I m *n(0, -oo) |0) . 
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The function J represents a transition amplitude connecting the vacuum and a one- 
particle state. The eikonal line simulates interactions with partons moving in different 
directions: the direction n'^ is arbitrary, but off the light-cone, in order to avoid spurious 
collinear singularities. 

The factorization formula (13) also involves the eikonal approximation to the partonic 
jet J, which we call the 'eikonal jet'. It is defined by 

(15) J (^^^^,as{fi% e) = (0| $0(00, 0) $„(0, -00) |0) , 



where the velocity vector Pi of each jet is related to the corresponding momentum pi by 
Pi = f^iQo/V^, with Qo a hard scale such that pi ■ Pj/Qo is of order one for all 

Both the partonic jet (14) and the eikonal jet (15) have soft divergences, as well as 
collinear divergences associated to their light-like leg; thus, they display double poles 
order by order in perturbation theory. The double poles are however the same, since in 
the soft region J' correctly approximates J: singular contributions to the two functions 
differ only by hard collinear radiation. 

The final ingredient in (13) is the soft matrix, which we define by taking the eikonal 
approximation for all gluon exchanges: since soft gluons do not resolve the details of 
the hard interaction nor the internal structure of the jets, they couple effectively to 
Wilson lines in the colour representations of the corresponding hard external partons. 
Such exchanges mix the colour components of the amplitude, forming a matrix in colour 
space. Choosing a basis of independent tensors cl in color space, we write 

n 

(16) E(c^){a.}'5^^ (A •/3„«.(M'),e) = E n [*ft(o«.0)"..'/.] |0) • 

L {rik} »=1 

Note that in eq. (16) we keep all Wilson lines strictly on the light-cone {pf — and thus 
f3f = 0). Therefore, the soft matrix Slk displays not only single poles corresponding 
to large-angle soft gluons, but also double poles associated with overlapping soft and 
collinear singularities. Recall that the jet functions Ji also include the regions of over- 
lapping soft and collinear singularities. It is for this reason that in the factorization 
formula, eq. (13), we have divided each partonic jet Jj jet by Ji, thus removing from 
Ji its eikonal part, which is already accounted for in Slk- One observes then that the 
ratios Ji/Ji are free of soft singularities: they contain only single collinear poles at each 
order in perturbation theory. Similarly, the 'reduced' soft matrix 

(17) Slk {Pij,as{fi^),€' 



2 '"«V^ 



where 



(IS) <*-^^ 



2(A ■ Hi)' 2{p, ■ n,Y 
2 



n 



3 



10 



E. GARDI and L. MAGNEA 



is free of coUinear poles, and contains only infrared singularities originating from soft 
gluon radiation at large angles with respect to all external legs. The kinematic depen- 
dence of <S on pij will be explained below. 

5. — The cusp anomaly and the jet functions 

The third and final ingredient necessary for deriving the constraints on the singu- 
larities of on-shell amplitudes is the cusp anomaly. In the following we will recall the 
definition of the cusp anomalous dimension and explain its role in governing the kine- 
matic dependence of Wilson-line operators. This will allow us first to understand the 
structure of the cikonal jet function J to all orders in perturbation theory, and eventually 
to constrain the soft function <S. 

To this end, let us recall some general properties of operators that are composed 
of semi-infinite Wilson linos. The first observation is that all radiative corrections to 
such operators vanish identically in dimensional regularization, since the corresponding 
integrals involve no scale. This trivial result however involves cancellations between ultra- 
violet and infrared singularities; therefore, upon renormalization, J becomes non-trivial: 
the contribution of each graph equals minus the corresponding ultraviolet counterterm. 
As a consequence, using a minimal subtraction scheme, the result for J (or for S) at 
each order in is a sum of poles in e, without any non- negative powers. 

Let us now briefiy recall the renormalization properties of Wilson loops with cusps 
[6-11,52-54]. Consider first an operator 

(19) W(7i2,a.(M'),e) = (0| $„,(«., 0) (0, -oo) |0) , 

involving two semi- infinite rays, both off the lightconc {n\, ^ 0), which join at the 
origin to form a cusp with (Minkowski) angle 712, where 

(20) cosh (712)= ^"X^ - 

The contour closes at infinity and it is smooth everywhere except at the origin. Ref. [53] 
has shown that the presence of the cusp along the contour introduces an ultraviolet 
singularity which can be removed by a multiplicative renormalization constant, implying 
that 

(21) ^jj^ = - reusp(7i2, a.(/^')) = - a [712 coth(7i2) - 1] + O (a^) , 

where Cj = Ga or Cf depending on the representation of the Wilson lines. Considering 
the limit where ni or n2 is near the lightcone one finds that 

and 

^ ' rfln/i 2 \ \ ^ y ) ■ 
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It is clear that in the strictly light-like limit, the r.h.s. of (23) and (21) become singular. 
This is a collinear singularity, appearing on top of the ultraviolet singularity already 
present in W owing to the cusp. Therefore, if we consider directly the renormalization 
of the analogue of W with one of the rays being strictly lightlike, say n2 = — precisely 
the case of J in eq. (15) — we expect a singular anomalous dimension. Indeed, in 
dimensional regularization the renormalization group equation of J takes the form [39] 

^ ' ^djl — ^ ,as[n),ej = -7J-J — ,as[n ),e 

= 2 " )) In J " 4 ^"^(^ ' ' 

where the third term is singular, 0{l/e). Having seen the origin of this singularity in (23), 
it is not surprising that the dependence of 7^- on the kinematic variable 2(/3 • n)^ jr? — 
the second term in (24) — is governed by the same anomalous dimension, 7^-'(as), that 
governs the singularity of 7j^ — the third term in (24). This relation between kinematic 
dependence and singular terms, which we have now observed in 7^-. , is a general property 
of this class of operators which will be essential for what follows. 

To understand it from a different angle, let us now have another look at eqs. (21) 
and (24) considering the symmetry property of the eikonal Feynman rules (10) under 
rescaling of the eikonal velocity vectors. Clearly eq. (21) is consistent with this symmetry: 
any function of 712 defined in (20) would be. In contrast, in the strictly lightlike case 
of (24) , there is no kinematic variable that could be consistent with this symmetry. J can 
only depend on 2(/? • rif' jr? , as indeed can be confirmed by an explicit calculation, and 
therefore it breaks the rescaling symmetry: it depends explicitly on the normalization 
of /3. Note that rescaling of the vector n^, which is not light-like, remains a symmetry. 

Solving (24) we obtain a closed form expression for the eikonal jet [39], in terms of 
anomalous dimensions which depend just on the coupling. 



Ji 



(25) - -7^ ^a,(A , e)j In 

We sec that the entire kinematic dependence of J is associated with the breaking of 
the rescaling symmetry with respect to the lightlike direction /3^; it is directly related 
to presence of double poles in J , and it is governed by the cusp anomalous dimension, 
l^K^OLs)- In the following, we will show how this observation, made in [39], allowed us to 
constrain the kinematic dependence of the soft function iS. 

Before turning to the soft function, let us quote the equivalent expression for the 
partonic jet J^, which will be of use in the following. The partonic jet has an infrared 
singularity structure similar to the eikonal jet, however it has a finite ultraviolet anoma- 
lous dimension, the one we have encountered in (5), 



" J j \ 2 9 



(26) Mxrl'^-^d 2 7 ,c.a(M^),e) ^-77.(a.(/x^)). 
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As a consequence, the partonic jet function can be written as 



(27) + ^ - -7if (a«(A , e)) In ^^^^ J + 2%("«(^ ' ^)) 

whore iJj is a finite coefficient function, independent of /i^. Eq. (27) displays the fact that 
in addition to the collinear singularities generated by jj. the jet function (f4) involves 
soft (eikonal) singularities; these are summarized by the second line of (27). Beyond 
the terms, one finds single pole terms governed by d!^K This function, which is 
defined in eqs. (4.7) and (4.9) in [39], does not depend on the spin of parton i and it 
has a maximally non-Abelian structure. For simplicity, we further assumed here that 
s!^^ admits Casimir scaling, d'^^ = Tj • Tj fc, although this may not hold beyond three 
loops (and it would not be important in what follows). In contrast, 7j. , which governs 
the collinear singularities, does depend on the spin of parton i (it differs for quarks and 
for gluons, see Appendix A of [41]) and it is not maximally non-Abelian. As anticipated, 
the double poles in (25) and (27) are the same, while single poles differ. 



6. — Derivation of the constraints on soft singularities 

We arc finally in a position to derive the promised constraints on soft singularities. We 
will show, in particular, that the relation established above, considering the case of the 
eikonal jet, between kinematic dependence of single pole terms and the cusp anomalous 
dimension, generalises to soft singularities in multi-leg amplitudes. These singularities 
are described by the function S in (13). S is defined by (16) and it obeys a matrix 
evolution equation of the form 

(28) {Pi ■Pj,a„e) = -Y^ Vjj (A • , a„ e) Sjk {Pi ■ Pj ,a„e) . 

The soft anomalous dimension matrix depends on all the kinematic invariants in 
the process, and it is a priori a very complicated object. It encapsulates the correlation 
between colour and kinematic degrees of freedom, which may be of increasing complexity 
as one considers higher loop corrections (fig. 2). 

We recall that in S all the Wilson lines are lightlike, /3| = 0, Vz. Therefore, similarly 
to J, we expect this function to break the rescaling invariance with respect to each of the 
velocities Pi. This was already taken into account in assigning the arguments in S and 
in F'^ : these functions depend on the set of Lorentz invariants Pi ■ Pj , and thus violate 
the rescaling symmetry. We will be able to constrain - and thus <S — because we 
know exactly how this symmetry is violated. 

The key point is that the amplitude M itself cannot depend on the normalization 
one chooses for the velocities appearing in eikonal functions. Thus, in the factorization 
formula, eq. (13), this dependence must cancel out. This cancellation can only involve 
the eikonal functions S and J7i, not the partonic jet or the hard function, which depend 
directly on the dimensionful kinematic variables Pi . The form of the factorization formula 
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(13) implies in fact that the cancellation of any rescaling violation must occur within the 
reduced soft function <S, defined in (17). This is intuitively clear: we saw that rescaling 
violation is intimately related to the presence of double poles, and that both are governed 
by the cusp anomalous dimension The soft function <S, much like the eikonal jets, is 
defined with light-like Wilson lines, thus including regions of overlapping ultraviolet and 
coUincar singularities, which arc the origin of double poles as well as rescaling violation 
at the single pole level. Upon dividing S by the product of all eikonal jets, as done in 
(17), these regions are removed, yielding S, which describes large-angle soft singularities, 
and is entirely free of double poles and of the associated violation of rescaling symmetry 
at the single pole level. Given the kinematic dependence of iS and Ji, and the expected 
recovery of the symmetry /?j — > Ki0i, Vi, we deduce that S can only depend on the 
variables pij, defined in (18). 

To proceed, it is useful to consider the renormalization group equation for the reduced 
soft function. In analogy with (28) we have 



(29) 



[pij, as, = - X] ^^'J' ^■^'^ (Pij^ois, e) 



In contrast to , the anomalous dimension of the reduced soft function, T^, is finite (<S 

itself has only single poles) and invariant under rescalings, as reflected in the fact that it 
must depend on Pi though pij only. Using the definition of the reduced soft function in 

eq. (17) we can directly relate the anomalous dimension to and to the anomalous 
dimension of the eikonal jets, of eq. (24). We obtain 



rfj {Pij, as) = Tjj (A • Pj,as, e) - Su ^ 



2(/3fc • Uk 
iJy I ^ — — ,a,,e 



(30) 



n r 



= rfj(A-/3,-,a«,e)-(5,jX 



k=l 



(2{pk-nkY\ 


1 /■ 


V nl ) 





^2 \ 



This equation implies highly non-trivial constraints on soft singularities. It tells us pre- 
cisely how the double poles and the rescaling violation of the single poles in cancel 
out. In particular, observing that the jet terms are diagonal in colour space (they 
are proportional to the identity matrix), we deduce that 

• off diagonal terms in must be finite, and must depend only on conformal cross 

ratios, 



(31) 



Pijkl 



(A-/3,)(/3fc.A) 

{Pi ■ PkXPj ■ Pi) 



Pij Pkl 
Pik Pjl 



1/2 



-i7r(Ai3 -I- Afci —^ik—^jl) 



which can be interchangeably expressed in terms of Pi ■ Pj (the arguments of F*^), 
or in terms of pij (the arguments of F"^); 
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• diagonal terms in have a singularity determined by jk, according to 



(32) 




1 




(a,(^2^e))+0(e°) 



and must contain finite terms depending on Pi ■ Pj in a way tailored to combine 
with the {Pi ■ UiY /n1 dependence of the various jet functions to generate pij. 

The constraints on the structure of the anomalous dimension can be compactly 
expressed by taking a logarithmic derivative of eq. (30) with respect to {Pi ■ /n1. On 
the l.h.s one uses the chain rule: for any function F which depends on {Pi ■ rii)^ jnl only 
through the combinations of (18), one has 



On the r.h.s of cq. (30), the derivative with respect to {Pi ■ nij^ /n1 acts only on the 
corresponding term. The resulting equations are 



Thus, there are n constraints for an n legged amplitude, each of which is a matrix 
equation (holding for each matrix element {I, J)). This set of constraints holds in any 
colour basis, and to all orders in perturbation theory. Its most intriguing aspect is that 
it correlates the kinematic dependence of the (reduced) soft matrix with its dependence 
on the colour degrees of freedom: the l.h.s in (34) is a siun of non-diagonal matrices in 
colour space, while the r.h.s is proportional to the identity matrix. 

7. — Solving the equations 

Given n independent equations and n{n — l)/2 kinematic variables it is clear at the 
outset that eq. (34) alone is not sufficient to uniquely fix the kinematic dependence of 

in the multi-leg case. For n = 2,3 eq. (34) does have a unique solution (sec sec. 4 
and Appendix A in [39]). This is already an important step, extending previously known 
results for the singularity structure to all loops. For n > 4 partons, however, the number 
of kinematic variables exceeds the number of equations, and additional constraints will be 
needed. Nevertheless, we will see that a minimal solution, consistent with all information 
known to date, naturally emerges out of eq. (34). 

Considering eq. (34), we note that depends implicitly on the colour representation 
of parton i. To solve the equations we need to make this dependence explicit. Given 
that 7^^' admits Casimir scaling (9) at least to three loops, we write 



(33) 




(34) 




Vi, /, J . 



(35) 



n^K = {as) + Ik 
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where = 0{ag) accounts for possible dependence on the representation of parton 
i through higher-order Casimir operators. It is presently an open question(^) whether 
such terms appear. 

Our constraints now take the form 



(36) 2^^Tr7:-Tr*(Pii,a.) 



Using the linearity of these equations we can obviously write the general solution as a 
superposition of two functions 

(37) r^ = r|c. + rl.c. 

which are, respectively, solutions of the equations 

(38) E ah|^ ^Q-C- ^P'J' = I • T, 7^ (a,) , Vi , 



(39) Edi^)^^-C.iPiJ,c.s) = \l'^^{as), Vi. 

Here Q.C. and H.C. stand for Quadratic Casimir and Higher-order Casimir, respectively. 
Let us now focus on determining Tq ^ , leaving aside T^ q = 0(af), which will be briefly 
discussed in sec. 8. 

A solution for , obeying eq. (38), is given by 

- 1 1 - " 

(40) {pij,as) = - g 7K (««) E 1^('^^^) • + 2 ^5("«) E • , 

where ^^^^ term in (40) indiactes a sum over all pairs of hard partons, 

forming a colour dipole; each dipole is counted twice in the sum. Note this term carries 
the entire dependence on kinematics, correlating it with the colour structure. In contrast, 
the second term is independent of kinematics and is proportional to the unit matrix in 
colour space. Analysis of the n = 2 case (the Sudakov form factor) allows to identify this 
function (see eqs. (4.7) and (4.9) in [39]) as the one appearing at single-pole level in the 
partonic jet, the last term in eq. (27). 

It is easy to verify that (40) satisfies (38): taking a derivative with respect to ln(py ), 
for specific partons i and j, isolates the color dipole Tj • Tj; summing over j for fixed i, 
and enforcing colour conservation, given by eq. (6), one recovers cq. (38). 



(^) An argument against Casimir scaling has been made [55,56], based on the dependence of 

7^' on the representation in the strong coupling limit at large Nc- The argument is based on a 
class of antisymmetric representations with k indices, where the ratio k/Nc is kept fixed when 
iVc — » 00; in this case the strong coupling limit of 7^^ does not admit Casimir scaling. 
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Integrating the renormalization group equation (29), with given by eq. (40), we 
obtain an expression for the reduced soft function. 



S {pij,as,e) = exp 



(41) 







1 


1-5 


L A2 


2 



i=l 



^^K{as{X\e))Y,ln{pij)Ti-Tj I. 



Substituting eq. (41) into the factorization formula, eq. (13), together with the corre- 
sponding expression for the partonic jet, eq. (27), we obtain a complete description of 
the singularity structure of the amplitude, cq. (5). Note that the Z factor and the hard 
amplitude Ti. in (4) arc separately independent of the auxiliary vectors n^, as they must 
be. In contrast, the various elements in the factorization formula (13) do depend on these 
vectors. The cancellation of this dependence is non-trivial: it is guaranteed by the fact 
that S admits the constraints of (34), and by the fact that the kinematic dependence 
of the singularities of the partonic jet function (27) matches the one of the eikonal jet, 
eq. (25). It is essential that all the single pole terms that carry n, dependence in the 
various functions arc governed by the cusp anomalous dimension ^/k alone. Indeed, to 
obtain eq. (5), we combine terms proportional to in the soft and jet functions. In 
doing so we use colour conservation, X^^^j Tj = — Tj, as well as the relation between the 
kinematic variables of the various functions, 
(42) 



In 



(2p^ ■ Uif 



In 



(2Pi -rijf 



In 



\ 



\ 



Note also that the poles associated with S-g{as{X'^,e)) cancel out between the soft and 
jet functions, given by eqs. (41) and (27), respectively. 

It is interesting to compare at this point our approach to that of Becher and Neubert 
in Ref. [41]. The final expression at the amplitude level, eq. (5), is the same. The set of 
constraints, eq. (48) in [41], is also equivalent. The underlying factorization scheme, and 
consequently the arguments leading to these constraints, are however somewhat different. 
In particular, ref. [41] does not define jet functions using auxiliary Wilson lines (n^ in 
our formulation); instead, it keeps track of the jets through their mass, taking pi slightly 
off the light cone, ^ 0. In their formulation, the equivalent of eq. (42) takes the form 
(eq. (43) in [41]) 

' ^ V ' ^ V ' 

J* Jj s 



which is again realised owing to the fact that in each function the corresponding logarithm 
is governed by the cusp anomalous dimension. Recall that in our derivation the argument 
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of the reduced soft function is dictated by rescaling invariance; in contrast, in [41] the 
argument of the soft function is essentially dictated by power counting, and it is not 
invariant with respect to rescaling. 

8. — Possible contributions beyond the sum-over-dipoles formula 

As already mentioned, the ansatz of (40) does not in general provide a unique solution 
of the available constraints. Thus eq. (5), while consistent with all existing calculations, 
may still be missing some large-angle soft singularities beyond a certain loop order. Such 
further singularities are however strongly constrained both in their functional form and in 
their color structure. It is worthwhile emphasizing that calculations in the large- A^c limit 
cannot resolve the question at hand, since planar eikonal diagrams necessarily factorize 
into a product of colour dipoles, those made of adjacent Wilson lines, and are therefore 
automatically consistent with the sum-over-dipoles formula. The analysis must therefore 
be done at finite Nc- Ref. [39] has identified two classes of corrections that may appear, 
and although some progress was made, neither of the two can be excluded to all orders 
at present. 

The first class of corrections corresponds to potential higher-order Casimir contri- 
butions. In case higher-order Casimir operators do show up in the cusp anomalous 
dimension at some loop order, i.e. jk in (35) does not vanish, the anomalous dimension 
of the reduced soft function of any amplitude will receive additional corrections. These 
corrections are subject to the very stringent constraints of (39). For amplitudes with two 
or three legs these corrections still have a dipolc structure (sec e.g. Appendix A in [39]), 
however, for amplitudes with four legs or more, non-trivial structures that couple more 
than two hard patrons may arise. 

The second class of corrections, which may be present even if "/k admits Casimir 
scaling, is given by solutions of the homogeneous equation associated with eq. (38). 
Indeed, adding to our ansatz any function A^{pij) satisfying 

(44) Eahfe^^(P^.'"«) = V^, 

one obtains a new solution of eq. (38). Eq. (44) is solved by any fimction of the confor- 
mal invariant cross ratios defined in (31). Any such solution has the property of being 
invariant with respect to velocity rescalings without involving the jets. Such functions 
can of course be written directly in terms of the original kinematic variables pj • pj and 
are therefore not constrained by soft-coUinear factorization. 

Interesting examples for A*^ in the four parton case were proposed in [39]: 

(45) E E i T,"T^Tf In (p,,fei) In (p^^y ) In {pa^H) , 

j,k,l a,b,c 

(46) E E "^-bc T]TlTf In^ (pijki) In' {Piklj) In' {pujk) , 

j,k,l a,b,c 

where the sum over partons is understood to exclude identical indices, and where colour 
conservation, Tf ~ — T^— T^— Tf , has been taken into account. Note that these functions 
are, by contruction, symmetric under the exchange of Wilson lines (Bose symmetry): this 
correlates colour and kinematic degrees of freedom. These functions, moreover, do not 
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contribute in the limit where any two hard partons become colhnear, and therefore they 
cannot be excluded using the properties of the splitting amplitude discussed in [41] . 

Functions of conformal-invariant cross ratios such as (45) correlate colour and kine- 
matic degrees of freedom of four partons. They cannot arise at two loops because two- 
loop webs can connect at most three partons. This explains, a posteriori, the findings 
of rcf. [38], which explicitly showed that there are no new correlations generated at the 
two loop order beyond those of pairwise inter actions (**). Non-trivial corrections to the 
the sum-over-dipoles formula can therefore first arise at three loops. 

Unfortunately, at three loops no complete calculation is available yet. However, sev- 
eral important steps have been taken. First, as already emphasized in [39], three-loop 
corrections to this formula must satisfy eq. (44) — they must be functions of confor- 
mal invariant cross ratios, and they must vanish identically in amplitudes of less than 
four legs. Beyond that, it was explicitly shown that the class of three-loop diagrams 
containing matter loops is consistent with the sum-over-dipoles formula [64]. 

A further step was taken in Rcf. [41], where it was shown that if dependence on the 
kinematic variables is assumed to be single-logarithmic, there is no possible structure 
that could appear at three-loops beyond the sum-over-dipoles formula. This argument is 
based on eliminating all possible structures using the factorization constraints discussed 
above, together with Bose symmetry, and an additional constraint on the singularity 
structure in the limit where two hard partons become coUinear based on the properties 
of the splitting amplitude. It should be emphasized that the assumption of single- 
logarithmic kinematic dependence is crucial for this argument, and this assumption may 
well be violated. It thus remains an open question whether new structures appear in the 
soft anomalous dimension at three loops. 

9. — Conclusions 

We have reviewed recent exciting progress in determining the infrared singularities of 
on-shell scattering amplitudes in massless non-Abelian gauge theories. It is now firmly 
established [39, 41] that the cusp anomalous dimension has a central role in governing 
soft singularities of multi-leg amplitudes with an arbitrary number of legs and for a 
general Nc. This role is summarized by a set of differential equations (34) constraining 
the kinematic dependence of the soft anomalous dimension matrix of any amplitude, to 
any loop order, in an arbitrary colour basis. These constraints are a direct consequence 
of factorization and of the special properties of soft gluon interactions with massless hard 
partons. 

The simplest solution to this set of constraints yields a closed form expression for the 
singularities of any massless scattering amplitude, eq. (5). According to this formula the 
correlations induced by soft gluon interactions between colour and kinematic degrees of 
freedom take the form of a sum over colour dipoles. No new correlations are generated 
by multi-loop webs (fig. 2): the colour matrix structure remains the same as at one loop, 
and the cusp anomalous dimension alone governs all non-collinear singularities. 

We have further shown that possible corrections to this simple sum-over-dipoles for- 
mula belong to one of two categories: ones that are generated by potential higher-order 
Casimir contributions to the cusp anomalous dimension, which must then satisfy eq. (39), 



Note that new structure does appear in the case of scattering involving heavy quarks, as 
shown in refs. [57-63]. 
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and ones that can be written in terms of conformal invariant cross ratios (31), solving 
the homogeneous equations (44). The former may contribute to any amphtude starting 
from four loops, while the latter can only appear in amplitudes with four or more hard 
partons, starting at three loops. So far all explicit calculations are consistent with the 
sum-over-dipoles formula, but it remains an open question whether such corrections do 
show up at some loop order. 
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